Phase diagram of underdoped cuprate superconductors: 
effect of Cooper-pair phase fluctuations 
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In underdoped cuprates fluctuations of the phase of the superconducting order parameter play a 
role due to the small superfluid density. We consider the effects of phase fluctuations assuming the ex- 
change of spin fluctuations to be the predominant pairing interaction. Spin fluctuations are treated 
in the fluctuation-exchange approximation, while phase fluctuations are included by Berezinskii- 
Kosterlitz-Thouless theory. We calculate the stiffness against phase fluctuations, n s (uj)/m* , as a 
function of doping, temperature, and frequency, taking its renormalization by phase fluctuations 
into account. The results are compared with recent measurements of the high-frequency conductiv- 
ity. Furthermore, we obtain the temperature T* , where the density of states at the Fermi energy 
starts to be suppressed, the temperature T* , where Cooper pairs form, and the superconduct- 
ing transition temperature T c , where their phase becomes coherent. We find a crossover from a 
phase-fluctuation-dominated regime with T c oc n s for underdoped cuprates to a BCS-like regime for 
overdoped materials. 

PACS numbers: 74.20.Mn, 74.40.+k, 74.72.-h 



I. INTRODUCTION 

For about fifteen years cuprate high-temperature su- 
perconductors (HTSC's) have stimulated significant ad- 
vances in the theory of highly correlated systems as well 
as in soft condensed matter theory. Nevertheless, we still 
do not fully understand the various phases of these mate- 
rials. Of particular interest is the underdoped regime of 
hole-doped cuprates, in which the hole density (doping) 
x in the CuC>2 planes is lower than required for the maxi- 
mum superconducting transition temperature T c . In this 
regime the superfluid density n s decreases with decreas- 
ing doping and is found to be proportional to T C B Above 
T c , one finds a strong suppression of the electronic den- 
sity of states close to the Fermi energy, i.e., a pseudogap, 
which appears to have the same symmetry as the super- 
conducting gaprl Furthermore, there may be fluctuating 
charge and spin modulations (stripes). u 

It has been recognized early on that the small su- 
perfluid density n s leads to a reduced stiffness against 
fluctuations of Jrhp. phase of the superconducting or- 
der parameter .DtrQ Phase fluctuations are addition- 
ally enhanced because they are canonically conjugate 
to charge density fluctuations, which are believed to 
be suppressed. Furthermore, the cuprates consist of 
weakly coupled two-dimensional (2D) CuC>2 planes so 
that fluctuations are enhanced by the reduced dimension- 
ality Phase fluctuations might destroy the long-range 
superconducting order, although there is still a conden- 
sate of preformed Cooper pairs. In conventional, bulk su- 
perconductors this mechanism is not relevant, since the 
large superfluid density leads to a typical energy scale of 
phase fluctuations much higher than the superconduct- 
ing energy gap A, which governs the thermal breaking 
of Cooper pairs. Thus in conventional superconductors 
the transition is due to the destruction of the Cooper 
pairs and T c is proportional to A.Q On the other hand, 



the observation that T c oc n s in underdoped cupratesB 
indicates that the phase fluctuations drive the transition 
in this regime. The Cooper pairs only break up at a 
crossover around T* > T c . If the feedback of phase fluc- 
tuations on the local formation of Cooper pairs is small, 
T* is approximately given by the transition temperature 
one would obtain without phase fluctuations. Between T c 
and T* Cooper fjaica. exist but the order parameter is not 
phase coherent .QErQa'El Recent thermal-expansion exper- 
iments strongly support this pictureJilj However, there 
is no close relation between our T* and the mean-field 
transition temperature of Ref. [l0| which is determined 
by extrapolation from the low-T behavior of the expan- 
sivity. 

There is a third temperature scale T* with T* > T*, 
below which a pseudogap starts to open up as seen in 
nuclear mag n e ti c reson ance, tunneling, and transport 
experimentsJliHlS'tSlESEj It seems unlikely that the pseu- 
dogap at these temperatures is due to local supercon- 
ductivity. Rather, it is thought to be caused by .spin. 
fluctuationsQ or the onset of stripe inhomogeneities.E£lEZl 
Recept experiments on the Hall effect in GdBa2Cu307_,$ 
filmsEj also support the existence of two crossover tem- 
peratures T* and T*. In this work we are mostly con- 
cerned with the strong pseudogap regime T c < T < T*. 

Due to the layered structure of the cuprates, they be- 
have like the 2D XY model except in a narrow criti- 
cal range around T c , where, they show three-dimensional 
(3D) XY critical behavior .t£lE£l The standard theory for 
the 2D XY model, the Berezinskii-Kosterhtz-Thouless 
(BKT) renormalization group theory,EjE3Ei20 should 
thus describe tftgss. materials outside of the narrow crit- 
ical range .oo'E3i£a'E3E3 Also, recent transport measure- 
ments for a gate-doped cuprateEa with only a single su- 
perconducting Cu02 plane show essentially the same 
doping dependence of T c as found for bulk materials. 
BKT theory predicts a transition at a temperature T c < 



2 



T*, due to the unbindung of fluctuating vortex- antivortex 
pairs in the superconducting order parameter. Gaussian 
phase fluctuations are less important, since they do not 
shift T c .Lj In addition, the coupling of the phase to the 
electromagnetic field causes them to be gapped at the 
plasma frequency,E3 whereas the BKT picture of vortex 
unbinding remains basically unchanged.^ 



In the early days of HTSC's, BKT theory was en- 
voked to Autf r.prjjft . r^h number of experiments on bulk 
samples. Hl3EZIC3cj'EJ Recently, two experiments have 
lent strong additional support to the BKT description: 
First, Corson et aln3 have measured the complex conduc- 
tivity of underdoped Bi2Sr2CaCu20s+5 and extracted 
the frequency-dependent phase stiffness from the data. 
The authors interpret thw data in terms of dynamical 
vortex-pair fluctuationsE3'c3 and conclude that vortices — 
and thus a local superconducting condensate — exist up to 



at least 100 K_We discuss this assertion in Sec. $1% Sec- 



ond, Xu et al£3 have found signs of vortices at tempera- 
tures much higher than T c in underdoped L^-zSr^CuO^ 
in measurements of the Nernst effect. A recent reanalysis 
of the dataE-2 yields an onset temperature of vortex effects 
of 40 K for an extremely underdoped sample (x = 0.05) 
and of 90 K for x = 0.07. 



So far, we have not said anything about the supercon- 
ducting pairing mechanism. There is increasing evidence 
that pairing is mainly due to the exchange of spin fluc- 
tuations. The|jpnsfl™iaa_fluctuation-exchange (FLEX) 
approximatioriJo'^-cZrEacj based on this mechanism de- 
scribes optimally doped and overdoped cuprates rather 
well. In particular, the correct doping dependence and 
order of magnitude of T c are obtained in this regime. On 
the other hand, the FLEX approximation does not in- 
clude phase fluctuations and we believe this to be the 
main reason why it fails to predict the downturn of T c in 
the underdoped regime. Instead, T c is found to approxi- 
mately saturate for small doping x. However, the FLEX 
approximation is able to reproduce two other salient fea- 
tures of underdoped cuprates, namely the decrease of n s 
and the opening of a weak pseudogap at T* , as we show 
below. 

This encourages us to apply the following description. 
We employ the FLEX approximation to obtain the dy- 
namical phase stiffness n s (u>) /m* , where n s (u>) is the gen- 
eralization of the superfluid density for finite frequencies. 
The static density n s (0) starts to deviate from zero at 
the temperature where Cooper pairs start to form and 
which we identify with T*. Then, phase fluctuations are 
incorporated by using the phase stiffness from FLEX as 
the input for BKT theory, which leads to a renormal- 
ized nf < n s and predicts a reduced T c . Then, we 
consider thq-dynamical case w > and use dynamical 
BKT theory23c3 to find the renormalized phase stiffness 
nf*(uj)/m* and compare the results with experimentsMJ 



II. STATIC CASE 

Transport measurements for a gate-doped cupratefEl 
show that the superconducting properties are determined 
by a single CuC-2 plane. The simplest model believed to 
contain the relevant steong correlations is the 2D one- 
band Hubbard model.E3 We here start from the Hamil- 
tonian 



H 



(4< 



C ja Ci, 



U y^n iT n a . (1) 



Here, c| CT creates an electron with spin a on site i, U 
denotes the on-site Coulomb interaction, and is the 
hopping integral. Within a conserving approximation, 
the one-electron self-energy is given by the functional 
derivative of a generating functional <&, which is related to 
the free energy, with respect to the,dressed one-electron 
Green function Q, E = S^[H]/SgiB On the other hand, 
the dressed Green function is given by the usual Dyson 
equation Q^ 1 = Q^ 1 — E in terms of the unperturbed 
Green function (Jo of the kinetic part of H alone—These 



equations determine the dressed Green functkt 
The T-matrixB or FLEX approximationlfcl 
is distinguished by the choice of a particular infinite sub- 
set of ladder and bubble diagrams for the generating 
functional The dressed Green functions are used to 
calculate the charge and spin susceptibilities. From these 
a Berk-Schrieffer-typd£3 pairing interaction is contructed, 
describing the exchange of charge and spin fluctuations. 
In a purely electronic pairing theory a self-consistent de- 
scription is required because the electrons do not only 
form Cooper pairs but also mediate the pairing inter- 
action. The quasiparticle self-energy components X v 
(v = 0, 3, 1) with respect_to the Pauli matrices r„ in the 
Nambu representation,0E3 i.e., Xq — oj(1 — Z) (renor- 
malization), X% — £ (energy shift), and X\ = <j> (gap 
parameter) are given by 

x u {\t,u) = — V / dfi [p fl (k-k',n)±p c (k-k',fi)] 

iV w Jo 



duj I(lu, f2, u> ) A„(k , oj ). 



(2) 



Here, the plus sign holds for Xq and X$ and the minus 
sign for X\. The kernel / and the spectral functions A v 
are given by 



A u (k,Uj) = 



f(-uj')+b(Q) f(u') + b(Q) 



u> + id — Q — 



0J' 



iS + Q-, 



1 a v (k,w) 
— 1m — 

7T -D(k,w) 



7.( 3 ) 
(4) 



where ao — ojZ, = £k + £, a\ = <j), and 
D= (uZ) 2 -[e k + tf-<j> 2 . 



(5) 

Here, / and b are the Fermi and Bose distribution func- 
tion, respectively. We use the bare tight-binding disper- 
sion relation for lattice constant a = b = 1 . 



£k = 2t (2 — cos k x — cos k y — fi) . 



(6) 
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The band filling n — l/^Ek^k i s determined 
with the help of the k-dependent occupation number 
rik = 2 duj f{ui) JV(k, w) which is calculated self- 
consistently. n — 1 corresponds to half filling. The inter- 
actions due to spin and charge fluctuations are given by 
P s = (27r)- 1 U 2 lm(3 Xs -Xso) ^ith X s = Xso (l~U X so) Jl 
and P c = (27r) _1 C/ 2 Im(3xc - Xco) with Xc = XcO (1 + 
Uxcq)^ 1 - In terms of spectral functions one has 



ImXsO,co(q, uj) 



N 



&J [f(oj')-f(oj' + u;)] 



[iV(k + q,u/ + w)JV(k,ci/) 

k 

±A 1 (k + q,uj' + Lu)A 1 (k,uj')]. 



(7) 



Here, N(k,w) — ^4o(k, uj) + Aa(\s., uj), and the real parts 
are calculated with the help of the Kramers-Kronig rela- 
tion. The substracted terms in P s and P c remove double 
counting that occurs in second order. The spin fluctua- 
tions are found to dominate the pairing interaction. The 
numerical calculations are performed on a square lattice 
with 256 x 256 points in the Brillouin zone and with 
200 points on the real uj axis up to 16 1 with an almost 
logarithmic mesh. The full momentum and frequency 
dependence of the quantities is kept. The convolutions 
in k space are carried out using fast Fourier transforma- 
tion. The superconducting state is found to have d x 2_ y 2- 
wave symmetry. T* is determined from the linearized 
gap equation. 

A field-theoretical derivation of the effective action of 
phase fluctuationa3'oE30 shows that the phase stiffness 
for frequency uj = is given by the 3D static superfluid 
density divided by the effective mass, n s (x,T)/m* . This 
quantity is given by 



with 



In,s — 



{In - Is) 



dio <Ji' S (uj), 



(8) 



(9) 



where af 1 (erf) is the real part of the conductivity in the 
normal (superconducting) state. Here we utilize the /- 
sum rule J Q dioa\(uj) — ire 2 n/2m* where n is the 3D 
electron density. The interpretation of Eq. (g) is that 
the spectral weight missing from the quasi-particle back- 
ground in u(uj) for T < T* must be in the superconduct- 
ing delta-function peak. 

<t(uj) is calculated in the normal and superconducting 
states using the Kubo formulaEao 

*(u>) = -=--/ <b>'[f(u/)-f (<*>'+«>)] 

he uj J_ 00 

k 

+ A 1 (k,uj'+uj)A 1 (k,uj')], (10) 
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FIG. 1: Static superfluid density as a function of temperature 
for three values of the doping x (symbols). The solid curves 
are fits of power laws with logarithmic corrections as explained 
in the text. The intersection of n a (T)/m* with the dashed 
line represents a simplified criterion for the BKT transition 
temperature T c . 



where t>k,i = de^/dki arc the band velocities within the 
CuC>2 plane and c is the c-axis lattice constant. Vertex 
corrections are neglected. 

The superfluid density (phase stiffness) n s /m* ob- 
tained in this way is shown in Fig. ^ for the three dop- 
ing values x — 0.091 (underdoped), x — 0.155 (approx- 
imately optimally doped), and x — 0.222 (overdoped). 
The figure also shows fits to the data at given dop- 
ing level, where we assume the form In n s (T) /to* = 
ao + Oi ln(T* — T) + a 2 In (T* — T) + . . ., i.e., a power-law 
dependence close to T* with logarithmic corrections. We 
use the fits to extrapolate to T = 0. The results show 
that T* depends on x only weakly in the underdoped 
regime but decreases rapidly in the overdoped. We come 
back to this below. Furthermore, n s /m* increases much 
more slowly below T* in the underdoped regime and ex- 
trapolates to a smaller value at T = 0. 

We have also calculated n s in units of the total hole 
density n, shown in Fig. ||, finding that n s /n is signif- 
icant reduced below unity, in agreement with experi- 
ments but in contradiction to BCS theory. The reduction 
is strongest for the underdoped case. Our results show 
that spin fluctuations can explain most of the observed 
reduction of n s . Also note that n s is linear in temper- 
ature for T — ► because of the nodes in the gap. The 
inset in Fig. | shows A 3 (0)/A 3 (T), where the penetra- 
tion depth isB A cx n s 1 , as a function of T* — T. The 
FLEX approximation yields A 3 (0)/A 3 (T) oc T* - T. The 
same power law has been found experimentally by Ka- 
mal et alt3 It has been attributed to critical fluctuations, 
starting about 10 K below the transition temperature,E2l 
since it coincides with the critical exponent expected for 
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FIG. 2: Ratio of superfluid density to total hole density for 
the same doping values x as in Fig. |l| The inset shows 



A 3 (T = 0)/A 3 (T) = nJ 2 (T)/n s r \T = Q), where A is the Lon 
don penetration depth, as a function of (T* — T). 



3/2, 



the 3D XY model. We here obtain the same power law 
from the FLEX approximation, which is purely 2D and 
does not contain critical fluctuations. Instead this rapid 
increase of n s oc 1/A 2 below T* is due to the selfcon- 
sistcncy, which leads to a more rapid opening of the gap 
than in BCS theory. We thus conclude that, while critical 
3D XY flu pt i a an ions are expected in a narrow tempera- 
ture range,tac2l they are not the origin of the observed 
power law on the scale of 10 K. 

Now we turn to the renormalization of n s due to phase 
(vortex) fluctuations. The BKT theory describes the un- 
binding-of thermally created pancake vortex-antivortex 
pairs£lrE3 The relevant parameters are the dimensionless 
stiffness K and the core-energy E c of vortices. The stiff- 
ness is related to n s byEJ 



K(T) = f3fi 



2 n s (T) d 



4' 



(11) 



where (3 is the inverse temperature and d is the aver- 
age spacing between Cu0 2 layers. Since we use a 2D 
model to describe double-layer cuprates, we set d to 
half the height of the unit cell of the typical represen- 
tative YBa2Cu306+ a . The stiffness K is also a mea- 
sure of the strength of the vortex-antivortex interaction 
V = iTikoTK ln(r/Vo). Here, tq is the minimum pair 
size, i.e., twice the vortex core radius, which is of the 
order of the in-plane Ginzburg-Landau coherence length 
£ a {,. For the core-energy we use an approximate result 
by Blatter et a/.,E3 E c — irksT K InK, where k is the 
Ginzburg parameter. Starting from the smallest vortex- 
antivortex pairs of size ro, the pairs are integrated out 
and their effect is incorporated b-^an approximate renor- 
malization of K and the fugacitylEl y = e~P Ea . This leads 



to the Kosterlitz recursion relations 



dl 

dK 

~dT 



-4ir 3 y 2 K 2 , 



(12) 
(13) 



where I = ln(r/ro) is a logarithmic length scale. For 
T > T c , K goes to zero for I — > oo, so that the interac- 
tion is screened at large distances and the largest vortex- 
antivortex pairs unbind. The unbound vortices destroy 
the superconducting, order and the Meifiner effect and 
lead to dissipation^ For T < T c , K approaches a finite 
value Kr = lim/^oo K and y vanishes in the limit I — > oo 
so that there are exponentially few large pairs and they 
still feel the logarithmic interaction. Bound pairs reduce 
K and thus n s , but do not destroy superconductivity. At 
T c , Kr jumps from a universal value of 2/tt to zero. The 
values of T c shown below are obtained by numerically 
integrating Eqs. Jl2] ) and (13) with n s taken from an in- 
terpolation between the points in Fig. [j]. It turns out 
that the renormalization of K for T < T c is very small 
so that one obtains T c from the simple criterion 



K{Tc) 



2 

7T 



n s (T c ) 2 4fc B T c 



7T h 2 d 



(14) 



for the unrenormalized stiffness with an error of less 
than 1 %. Eq. (|lj) is satisfied at the intersection of the 
n s (T)/m* curves with the dashed straight line in Fig. [I]. 
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FIG. 3: Temperature scales of the cuprates as functions of 
doping x. T c (solid circles) is the transition temperature 
obtained from the FLEX approximation with phase fluctua- 
tions included by means of BKT theory. At T* (open circles) 
Cooper pairs start to form locally; this temperature is given 
by the transition temperature obtained from the FLEX ap- 
proximation with spin fluctuations alone. The crosses show 
the superfluid density (phase stiffness) n s (T = 0) /m* for com- 
parison. This curve has been scaled so that it agrees with T c 
in the underdoped regime. 
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From BKT theory we obtain two important quantities: 
the transition temperature T c and the renormalized stiff- 
ness Kr, which determines the renormalized superfiuid 
density (phase stiffness) 



300 



f3h 2 d 



Kr- 



(15) 



In Fig. |^ we plot the transition temperature T c and the 
temperature T* where Cooper-pairs form. For decreas- 
ing doping x, T* becomes nearly constant and even de- 
creases for the lowest doping level, consistent with the 
strong decrease of the onset temperature of v pr . tex ef- 
fects at even lower doping found by Xu et aZ.HllJ We 
have also calculated der superconducting gap A extra- 
polated to T = (not shown). Ao is here defined as half 
the peak-to-peak separation in the density of states. We 
find approximately Ao oc T*. 

Phase fluctuations lead to a downturn of T c in the un- 
derdoped regime. However, this reduction is not as large 
as experimentally observed and our value x ~ 0.14 for 
the optimal doping is accordingly smaller than the ex- 
perimental one ofiK 0.16.E3 We suggest that one origin 
of this discrepancy is the neglect of the feedback of phase 
fluctuations on the electronic properties. 

Figure ||| also shows the superfiuid density n s (0) /to* 
extrapolated to T = 0, scaled such that it approaches T c 
in the underdoped regime. The density increases approx- 
imately linearly with doping except for the most over- 
doped point, where it turns down again. This behavior 
agrees well with angl£-resolved photoemission (ARPES) 
results of Feng et aL 3 and with recent /iSK experiments 
of Bernhard et al£3 In Ref. [36] a maximum in n s at a 
unique doping value of x mSLX « 0.19 is found for vari- 
ous cuprates, while we obtain x max w 0r?0. Our results 
are consistent with the Uemura scalingcl T c oc n s (0) in 
the heavily underdoped regime and with the BCS-like 
behavior T c rs T* oc Ao in the overdoped limit. T c in- 
terpolates smoothly between the extreme cases. We find 
T c < T* even for high doping, since n s (T) and K(T) con- 
tinuously go to zero at T* so that Eq. ( jl4| ) is only satisfied 
at a temperature T c < T* . The results for the overdoped 
case may be changed if amplitude fluctuations of the cp. 
der parameter and their mixing with phase fluctuationsEll 
are taken into account. Amplitude fluctuations are gov- 
erned by A, which becomes smaller than the energy scale 
of phase fluctuations in the overdoped regime. 

The situation is complicated by the Josephson cou- 
pling between Cu02 layers. This coupling leads to the 
appearance of Josephson vortex lines connecting the pan- 
cake vortices between the layers. E2I They induce a lin- 
ear component in the vortex-antivortex interaction. This 
contribution becomes relevant at separations larger than. 
A = d/e, where e < 1 is the anisotropy parameter .E3 
A acts as a cutoff for the Kosterlitz recursion relations 
and eventually leads to an increase of T c relative to the 
BKT result T C BKT and-,te,|%j-|breakdown of 2D theory 
close to the transition.Effial2§EZl The experiments of Cor- 
son et clIeA also show that the BKT temperature T® KT 
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FIG. 4: Temperature T* at which a small suppression of the 
density of states at the Fermi energy (weak pseudogap) ap- 
pears. The temperatures T£ and T c from Fig. ^ are also 
shown. The inset shows the suppression of the density of 
states (in arbitrary units) for x = 0.155 and T — 4.5 T* (solid 
line), T = 2.3 T c * w T* (dashed line), and T = 1.01 T c * (dotted 
line) . 



extracted from the data is significantly smaller than the 
experimental T c . Thus T c as calculated here is a lower 
bound of the true transition temperature. 

The feedback of phase fluctuations on the electrons is 
not included in our approach. We expect the phase fluc- 
tuations in this regime to lead to pair breaking How- 
ever, simulations of the XY model suggest that this feed- 
back is rather weak.Ea Neglecting the feedback, the elec- 
tronic spectral function shows the unrenormalized super- 
conducting gap for T c < T < T*. Since there is no su- 
perconducting order in this regime, we identify this gap 
with the (strong) pseudogap, which thus is automatically 
d 2 .2„j / 2-wave-like and of the same magnitude as the su- 
perconducting gap for T < T c . Thus in this picture the 
pseudogap is due to local Cooper pair formation in the 
absence of long-range phase coherence. Pair breaking due 
to phase fluctuations should partly fill in this gap. 

Figure [| shows T c , T*, and T* on a different tem- 
perature scale. T* is the highest temperature where a 
weak pseudogap is obtained from FLEX, i.e., where the 
density of states at the Fermi energy starts to be sup- 
pressed. The inset shows this suppression for x = 0.155. 
The temperature T* becomes much larger than T c in-Lhc 
underdoped regime, in agreement with experiments.Ej 

To conclude this section, we discuss the effect of a 
normal-state pseudogap due to a mechanism other than 
incoherent Cooper pairing. Let us assume a suppres- 
sion of the density of states close to the Fermi surface in 
the normal state, e.g., due to the formation of a charge- 
density wave.O This decreases the number of holes avail- 
able for pairing and should thus reduce T c . To check 
this, we have performed FLEX calculations with a pseu- 
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FIG. 5: Transition temperatures in the presence of a normal- 
state pseudogap. The open squares show the transition tem- 
perature T* obtained from the FLEX approximation with a 
d-wave pseudogap in the normal-state dispersion. rBftei am- 
plitude of the pseudogap is taken from experiments. Ej'Cj The 
open circles show the corresponding values without a pseudo- 
gap, see Fig. [| The solid squares denote T c in the presence 
of the pseudogap and with phase fluctuations included, as- 
suming the two effects to be independent. The solid circles 
show the corresponding results without pseudogap. The in- 
set gives the phase stiffness n 3 /m* for the doping x = 0.122 
with (lower curve) and without (upper curve) the pseudogap. 
Intersections with the dashed line give the simple criterion 
([ti]) for T c . One clearly sees that a normal-state pseudogap 
increases the effect of phase fluctuations due to the slow in- 
crease of n s /m* below T* . 



dogap of the form Ak = Ao^osfcr — cosk y ) included 
in the normal-state dispersion. The doping-dependent 
amplitude Ao is chosen in accordance with ARPES ex- 
periments by Marshall et alEB and by Ding et al\3 The 
results are shown by the open squares in Fig. The 
curve merges with the T* curve without pseudogap (open 
circles) at x = 0.155, since-kese the pseudogap is exper- 
imentally found to vanishOEj It is apparent that T* is 
indeed strongly reduced in the underdoped regime. Thus 
this density-of-states effect is a possible alternative ex- 
planation for the observed downturn of T c . 

Next, we consider phase fluctuations in the presence 
of a normal-state pseudogap. The T c values naively ob- 
tained from BKT theory for this case are shown in Fig. || 
as the solid squares. Phase fluctuations reduce T c even 
more, in particular for x = 0.122. This is due to the 
fact that the phase stiffness n s /m* increases much more 
slowly below T* in the presence of a pseudogap, as shown 
in the inset of Fig. ||, even if T* is only slightly reduced. 
The small stiffness makes phase fluctuations more effec- 
tive. However, in this picture the reduction of T c is prob- 
ably overestimated: Above, we have explained the pseu- 
dogap as resulting from incoherent Cooper pairing. This 



contribution to the pseudogap must not be incorporated 
into the normal-state dispersion to avoid double count- 
ing. This would increase the result for T c . It is clearly 
important to develop a theory that incorporates phase 
fluctuations, spin fluctuations, and possibly the charge- 
density wave on the same microscopic level. However, 
the inclusion of vortex fluctuations in a FLEX-type the- 
ory on equal footing with spin fluctuations would be a 
formidable taskJj 



III. DYNAMICAL CASE 

In this section, we calculate the dynamical phase stifL 
ness, which is the quantity obtained by Corson et alJt3 
We first note that the superfluid density can also be ob- 
tained from the imaginary part of the conductivity, 

^ = \ lim u<r§(w), (16) 

as can be shown with the help of Kramers-Kronig rela- 
tions. We have recalculated n s /m* in this way and find 
identical results compared to Eq. (||). 

The phase stiffness has also been obtainedpaL-oanzero 
frequencies using field-theoretical methods.SH'ES'LZl For 
small wave vector q — > 0, 



The imaginary part <rf (oj) of the dynamical conductivity 
is obtained from the FLEX approximation for the dynam- 
ical curr.eiLt-current correlation function using the Kubo 
formula.E3 For lo > one should not interpret n s (uj) as 
a density. Note also that n s (ui) is no longer propor- 
tional to the penetration depth of a magnetic field — for 
lj > there is also a contribution from the real part of 
the conductivity, i.e., the normal skin effect. 

The resulting phase stiffness n s {ui)/m* is shown in 
Fig. U for x = 0.122 (underdoped) at various temper- 
atures. At higher doping the results (not shown) are 
similar, only the typical frequency scale, which turns out 
to be the low-temperature superconducting gap Ao, is re- 
duced. We find a finite phase stiffness at ui > even for 
T > T*. At first glance this is surprising, since the phase 
is not well-defined for A = 0. Indeed, using a Ward iden- 
tity one can show that thq-Gaussian part of the phase 
action vanishes for T > T*E3 However, the phase action 
contains a contribution from the time derivative of the 
phase besides the stiffness term. While the total action 
vanishes, each term on its own does not. Thus the stiff- 
ness is finite but has no physical significance for T > T* . 

Even slightly below T*, n s (ui = 0)/m* obtains a sig- 
nificant finite value, leading to the Meifiner effect, and 
there is a considerable redistribution of weight from ener- 
gies roughly above twice the low-temperature maximum 
gap, 2Ao, to energies below 2Ao. This redistribution in- 
creases with decreasing temperature. Also, a peak devel- 
ops slightly below Aq followed by a dip around 2Aq, this 
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FIG. 6: Frequency-dependent phase stiffness n s (uj)/m* for 
doping x = 0.122 (underdoped) and temperatures ksT/t = 
0.012, 0.015, 0.016, 0.017, 0.018, 0.019, 0.0195, 0.02, 0.0205, 
0.021, 0.0215, 0.022, 0.0225, 0.023 [with decreasing n s (0)/m*]. 
t = 250 meV is the hopping integral. The frequency is given 
in units of t {% = 1). At T* w 0.023 t/k B = 66.5 K Cooper 
pairs start to form. Below T* there is a marked transfer of 
weight from energies above 2Ao to energies below, where Ao 
is the maximum gap at low temperatures as obtained from 
the FLEX approximation. 



Only vortex-antivortex pairs of size r < r w contribute 
to the screening. Hence, we cut off the renormalization 
flows at this length scale. To avoid an unphysical kink in 
nf(bj)/m* we use the smooth cutoff r 2 = r 2 + r§. 

The diffusion constant of vortices is not easy to calcu- 
late accurately. In the-absence of pinning, the theory of 
Bardeen and StepherJI3 yields 



D, 



27TC^>„ k B T 



bid 



(19) 



where c is the speed of light, £_ a b ~ r /2 is the coher- 
ence length, p n is the normal-state resistivity, (f>o = hc/2e 
is the superconducting flux quantum, and d is an effec- 
tive layer thickness. In the renormalization the quan- 
tity D®/tq enters, which according to Eq. fll9| ) is lin- 
ear in temperature. In the presence of a high density a£ 
weak pinning centers the diffusion constant becomesE-il 
D v = D° v exp(— E p /k B T), where E p is the pinning 
energy. Matters are complicated by the obsepwation 
that E p depends on temperature. Rogers et alEB find 
E p {T) fa E$(l- T/T*) with E°/k B « 1200 K for 
Bi2Sr2CaCu20g+5. Absorbing the constant term in the 
exponent into the prefactor, the result for the diffusion 
constant in natural units is 



Dr 



C v — — exp 



k B T 



(20) 



structure being most pronounced in the underdoped case. 
Since Ao is smaller in the overdoped regime, n s (Lo)/m* 
changes more rapidly for small u> in this case. It is of 
course not surprising that 2A is the characteristic fre- 
quency of changes in n s {oj)/m* related to the formation 
of Cooper pairs. 

We now turn to the question of how phase fluctuations 
affect the dynamical phase stiffness n s (u)/m*. This re- 
quires a dynamical generalization of BKTJ-lppry, which 
was first developed by Ambegaokar et aZ.oO Here, we 
start from a heuristic argumciiLfor the dynamical screen- 
ing of the vortex interactions An applied electromag- 
netic field exerts a force on the vortices mainly by induc- 
ing a superflow, which leads to a Lorentz force on the 
flux carried by the vortices. On the other hand, moving 
a vortex leads to dissipation in its core and thus to a fi- 
nite diffusion constant D v ,13 which impedes its motion. If 
one assumes a rotating field of frequency w, small vortex- 
antivortex pairs will rotate to stay aligned with the field. 
Large pairs, on the other hand, will not be able to follow 
the rotation and thus become ineffective for the screen- 
ing. A pair can follow the field if its component vortex 
and antivortex can move a distance 27rr during one period 
T^, = 2tt/lu. During this time a vortex can move a dis- 
tance of about the diffusion length \/D v T u — y / 2nD v /uj, 
so that the critical scale for the pair size is 



Dv 
2t:lu 



(18) 



where C v is a dimensionless constant. However, such a 
large value of Ep would lead to a sharp, step-like depen- 
dence of n s (<jj)/m* on temperature, in contradiction to 
the smooth behavior shown in Fig. 4 of Rcf . fill In view of 
these difficulties we treat D v /tq as a constant parameter 
and discuss the dependence on D v below. 

To find the effect of phase (vortex) fluctuations on the 
phase stiffness, the recursion relations (12) and dl3| ) are 
now integrated numerically up to the cutoff I = ln(r/ro), 
which depends on D v /tq. The resulting renormalized 
phase stiffness nf-(co)/m* for constant D v /vq = 10 17 s~ 1 
and x = 0.122 is plotted in Fig. [?]. Other values of D v 
give similar results. Of course, faster vortex diffusion 
shifts the features at given temperature to higher fre- 
quencies. The dashed lines denote the unrenormalized 
stiffness, i.e., the same data as in Fig. ^, albeit on an 
expanded frequency scale. The highest frequency used in 
Ref. |4| (600 GHz) corresponds to uj/t « 0.01, also indi- 
cated in Fig. 0. 

For T < T c (the upper six curves) the static renor- 
malization has been found to be small, see Sec. If. The 
renormalization at finite lo is even weaker so that the 
renormalized stiffness is in practice identical to the un- 
renormalized one, which has only a weak frequency de- 
pendence for low oj, in agreement with Ref. |4l] . 

When T is increased above T c (the lower five curves in 
Fig. 0), a strong renormalization of the stiffness due to 
phase fluctuations sets in starting at very low frequencies. 
The Meifiner effect is thus destroyed for all T > T c by 
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FIG. 7: Phase stiffness nf(Lo)/m* renormalized by vortex 
fluctuations for x = 0.122 at temperatures faT/t = 0.016, 
0.017, 0.018, 0.019, 0.0195, 0.02, 0.0205, 0.021, 0.0215, 0.022, 
0.0225 (heavy solid lines). The vortex diffusion constant has 
been chosen as D v /rQ = 10 17 s _1 . The unrenormalized stiff- 
ness is shown as dashed lines; these are the same data as in 
Fig. ^| Note the expanded frepliency scale. The highest fre- 
quency used by Corson et als3 is indicated by the vertical 
dotted line. 



FIG. 8: Renormalized phase stiffness nf(uj)/rn* for x = 0.122 
as a function of temperature for frequencies / = 100 GHz, 
200 GHz, 600 GHz (heavy solid lines). The unrenormalized 
stiffness is shown as dashed lines. The dotted line represents 
the approximate criterion Eq. (jwj) for the (u> = 0) BKT tran- 
sition. 



the comparatively slow vortex diffusion. With increas- 
ing temperature the onset of renormalization shifts to 
higher frequencies. At frequencies above this onset, the 
vortices cannot follow the field and thus do not affect the 
response, as discussed above. The onset frequencies are 
always much smaller than 2Ao. The features at the en- 
ergy scale 2Ao shown in Fig. || which are due to Cooper- 
pair formation, are unaffected by phase fluctuations and 
show no anomality at T c . They vanish only at T*. 

Finally, in Fig. ^| we plot the renormalized n^(u))/m* 
for x — 0.122 as a function of temperature for various fre- 
quencies. This graph should be compared to Figs. 2 and 
4 of Ref. |4l] — note that the quantity Tg given there is 
proportional to nf/m*. We note that Corson et alxA 
assume a thermally activated density of free vortices, 
nf oc exp(— Ef/T), for T not too close to the BKT tran- 
sition temperature, and a temperature-independent dif- 
fusion constantJZj Here, we instead integrate the recur- 
sion relations (|l2|) and (|l^) explicitly up to the dynamical 
length scale r so that we do not have to make an assump- 
tion onnj. One sees that even at / = 600 GHz the broad- 
ened BKT transition is still much narrower than found 
by Corson et alc3 From Eqs. (O) and ( p(i| ) it is clear- 
that the diffusion constant D v /tq increases with tem- 
perature. In the presence of pinning it increases rapidly 
around ksT ~ E®. Since a larger diffusion constant, i.e., 
more mobile vortices, leads to stronger renormalization, 
the transition in Fig. ^ would become even sharper if 
D v /r$ were an increasing function of temperature. 

Our results show that dynamical BKT theory together 



with Bardeen-Stephen theory for vortex diffusion and 
natural assumptions on pinning does not agree quan- 
titatively with the experimental results.EJ We conclude 
that the finite size effect apparent in the experimental 
data is not only due to the finite diffusion length. An- 
other possible source is the interlayer Josephson coupling, 
which leads to the apperance of the Josephson length A 
as an additional length scale, as discussed above.El This 
length scale leads to a cutoff of the recursion relations 
at I ~ ln(A/r ), which becomes small close to T* due 
to the divergence of ro ~ £ a 6 (neglecting the feedback of 
phase fluctuations on the quasiparticles). This broadens 
the transition but cannot easily explain the observed fre- 
quency dependence. On the other hand, the experimcp^ 
tal observation that the curves for various frequenciescll 
start to coincide where the phase stiffness agrees with the 
universal jump criterion (|l4|) supports an interpretation 
in terms of vortex fluctuations. We suggest that a bet- 
ter description of the interplay of vortex dynamics and 
interlayer coupling is required to understand the data. 

Note, the origin of the discrepancy may also lie in the 
FLEX results for n s (uj)/m*, which do not include all ef- 
fects of temperature-dependent scattering on the conduc- 
tivity a ,LJ and in the omission of the feedback of phase 
fluctuations on the electronic properties. Another effect 
neglected here is the possible coupling to a charge-density 
wave perhaps taking the form of dynamical stripes. 
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IV. SUMMARY AND CONCLUSIONS 

In the present paper we have obtained the character- 
istic energy scales of hole-doped cuprate superconduc- 
tors from a theory that includes both spin and Cooper- 
pair phase fluctuations. The former are described by the 
FLEX approximation, whereas the latter are included by 
means of the Bcrezinskii-Kosterlitz-Thouless (BKT) the- 
ory, taking the FLEX results as input. Phase fluctua- 
tions mainly take the form of vortex fluctuations, since 
Gaussian phase flucuations have a large energy gap. Vor- 
tices lead to the renormalization of the phase stiffness 
n s {uj)/m* to nf(u>)/m*. The stiffness at T — > shows 
a maximum at a dopintJevel of x ~ 0.2, in good agree- 
ment with experiments.Ej At the transition temperature 
T c the renormalized static phase stiffness nf(ui — 0)/m* 
vanishes, leading to the disappearance of the Meifiner ef- 
fect. The ideal conductivity is also destroyed by free vor- 
tices. T c is significantly reduced compared to the tran- 
sition temperature T* that would result from spin fluc- 
tuations alone. The T c determined from spin and phase 
(vortex) fluctuations is in much better agreement with 
experiments in the underdoped regime and shows a max- 
imum at optimum doping. Still, our approach does not 
explain the full reduction of T c . We believe that a further 
reduction of T c results from (a) the breaking of Cooper 
pairs by scattering with phase fluctuations and (b) other 
instabilities that reduce the density of states in the nor- 
mal state, for example a charge-density wave. Since the 
latter effect also suppresses n s /m*, phase fluctuations 



can become even more effective and reduce T c further. It 
would be desirable to include the pair-breaking effect of 
phase fluctuations and the possible formation of a charge- 
density wave, on the same microscopic level as the spin 
fluctuations.B 

For T c < T < T*, where phase-coherent superconduc- 
tivity is absent, phase fluctuations lead to a strong renor- 
malization of n s /m* at frequencies much smaller than 
2Ao. Our results show tie same trends as found in con- 
ductivity measurements^!! However, a three-dimensional 
description of vortex dynamics might be required to ob- 
tain a more quantitative agreement. Local formation of 
Cooper pairs still takes place in this regime. This leads to 
a strong pseudogap of the same magnitude Ao and sym- 
metry as the superconducting gap below T c . We also find 
a frequency dependence of n^(u)/m* at higher frequen- 
cies, u) > Ao, that is very similar to the superconducting 
phase. These features vanish only around T* . Finally, for 
T* < T < T* there is a weak suppression in the density of 
states at the Fermi energy. Our results reproduce several 
of the main features common to all hole-doped cuprate 
superconductors. We conclude that the exchange of spin 
fluctuations, modified by strong superconducting phase 
(vortex) fluctuations in the underdoped regime, is the 
main mechanism of superconductivity in cuprates. 
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